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ABSTRACT 


The  Navy's  Impact  Burial  Prediction  Model  (IMPACT  25)  determines  the  amount 
of  burial  a  mine  experiences  upon  impacting  the  marine  sediment.  Impact  burial 
calculations  are  derived  primarily  from  the  sediment  characteristics  and  from  the  mine's 
two-dimensional  air  and  water  phase  trajectories.  Accurate  burial  prediction  requires  that 
the  model's  air  and  water  phase  trajectories  reasonably  mimic  the  objects  tme  trajectory. 
IMPACT  25  assumes  that  the  objects  are  cylindrical  in  shape  and  calculates  the  air  and 
water  phase  trajectories  entirely  from  momentum  equations. 

In  order  to  determine  what  effect  a  varying  center  of  mass  has  on  a  mine's  water 
phase  trajectory,  a  Mine  Drop  Experiment  was  conducted.  The  experiment  consisted  of 
dropping  three  cylinders  of  various  lengths  into  a  pool  where  the  trajectories  were  filmed 
from  two  angles.  The  controlled  parameters  were,  the  ratio  of  mine  length  to  diameter, 
initial  velocity,  center  of  mass  position  and  drop  angle.  Results  indicate  that  center  of 
mass  position  has  the  largest  influence  on  the  object's  trajectory  and  that  accurate 
trajectory  modeling  requires  the  inclusion  of  both  momentum  and  moment  equations. 
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I.  INTRODUCTION 


On  December  31,  1991  the  Union  of  Soviet  Socialist  Republics  (USSR) 
effectively  ceased  to  exist  under  international  law  and  the  Cold  War  ended  (Fischer 
1999).  In  response,  the  Navy-Marine  Corps  team  developed  a  new  strategic  concept, 
“...From  the  Sea”  (FTS),  which  provides  a  framework  for  Naval  operations  into  the  21st 
century.  FTS  effectively  shifted  operational  focus  from  blue  water  operations  to  sea- 
based  power  projection  into  regional  littoral  areas.  FTS,  the  1994  revision  “Forward  ... 
From  the  Sea”  (FFTS)  and  “Operational  Maneuver  from  the  Sea”  (OMFTS)  all  provide 
guiding  principles  for  sea-based  power  projection  to  regional  littoral  areas  of  the  world. 

One  of  the  greatest  threats  to  U.S.  sea-based  power  projection  in  littoral  areas  is 
the  naval  mine.  Mines  were  first  developed  in  1776  and  have  been  used  in  most  major 
conflicts  since.  Today,  an  estimated  50  countries  possess  some  sort  of  mining  capability. 
(Lehr  2000)  Mines  can  be  used  in  both  offensive  and  defensive  roles.  Offensively,  they 
can  be  placed  in  enemy  waters  or  nearby  sea-lanes  in  order  to  harass  military  and 
commercial  shipping.  Defensively,  they  can  be  used  to  delay  or  prevent  amphibious 
assaults  or  to  deny  command  of  the  sea.  The  Wonsan  Korea  Mne  Crisis  and  Iraq’s  use  of 
mines  during  Desert  Storm  provide  excellent  examples  of  the  value  of  the  naval  mine  as  a 
defensive  weapon.  Shortly  after  the  October  1950  Wonsan,  Korea  mine  crisis,  then  Chief 
of  Naval  Operations  Admiral  Forest  Sherman  exclaimed, 

”...  when  you  can’t  go  where  you  want  to,  when  you  want  to,  you  haven’t 
got  command  of  the  sea.  Command  of  the  sea  is  the  bedrock  for  all  of  our 
war  plans.  We  have  always  been  submarine -conscious  and  air- conscious. 

We  have  now  commenced  to  become  mine- conscious...  beginning  last 
week.”  (Boorda  1999). 

Within  the  past  15  years  three  U.S.  ships,  the  USS  Samuel  B.  Roberts  (FFG-58), 
Tripoli  (LPH-10)  and  Princeton  (CG-59)  have  fallen  victim  to  mines.  Total  ship  damages 
were  $125  million  while  the  mines  cost  approximately  $30  thousand.  (Boorda  1999) 
Mines  have  evolved  over  the  years  from  the  dumb  “homed”  contact  mines  that  damaged 
the  Tripoli  and  Roberts  to  ones  that  are  relatively  sophisticated  -  non- magnetic  materials, 
irregular  shapes,  anechoic  coatings,  multiple  sensors  and  ship  count  routines.  Despite 
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their  increased  sophistication,  mines  remain  inexpensive  and  are  relatively  easy  to 
manufacture,  upkeep  and  place.  As  such,  they  are  an  efficient,  yet  potent,  force  multiplier 
and  are  widely  available  to  any  country  or  group  who  has  a  modest  ability  to  purchase 
them. 

Naval  mines  are  characterized  by  three  factors:  position  in  water  (bottom, 
moored,  rising,  floating),  method  of  delivery  (aircraft,  surface,  subsurface)  and  method  of 
actuation  (acoustic  and/or  magnetic  influence,  pressure,  contact,  controlled).  The  littoral 
battlespace  is  divided  into  five  regions  based  upon  water  depth.  Within  each  of  these 
regions  naval  forces  can  encounter  multiple  types  of  threats  (Fig.  1).  The  littoral  regions 
are: 

•  Deep  Water  (DW).  Water  depths:  >300  ft.  Threat:  mainly  moored  and 
rising  mines,  although  a  few  large  bottom  mines  exist. 

•  Shallow  Water  (SW).  Water  depths:  from  40  to  300  ft.  Threat:  bottom, 
moored  and  rising. 

•  Very  Shallow  Water  (VSW).  Water  depths:  from  10  to  40  ft.  Threat: 
bottom,  moored,  rising  and  controlled. 

•  Surf  Zone  (SZ).  Water  depths:  <  10  ft.  to  the  beach  itself.  Threat:  same  as 
VSW  but  land  mines  and  obstacles  can  also  be  encountered. 

•  Craft  Landing  Zone  (CLZ).  Water  depths:  the  beach  itself.  Threat: 
conventional  land  mines  and  obstacles.  (U.S.  Naval  Mine  Warfare  Plan 
2000) 

Hie  shift  in  focus  from  the  blue  water  to  the  littoral  has  brought  many  new 
challenges  to  the  warfighter.  The  greatest  is  what  impact  will  the  highly  variable  littoral 
environment  have  on  future  operations,  particularly  mine  countermeasures  (MCM).  The 
most  influential  environmental  parameter  to  successful  MCM  operations  is  the  local 
bathymetry  character  of  the  bottom.  This  key  parameter  often  determines  whether  an  area 
should  be  swept  or  hunted.  Bottom  clutter  in  the  form  of  rock  outcrops,  coral  reefs,  man¬ 
made  debris  and  irregularities  in  slope  provide  false  contacts  or  create  shadow  zones  that 
increase  overall  clearance  times.  Soft  bottom  sediments  such  as  marine  clays  and  silts 
cause  a  high  degree  of  mine  burial  upon  impact.  These  buried  or  partially  buried  bottom 
mines  are  of  greatest  concern  to  the  MCM  planner. 
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SEA  MINES  : 

•  Floating  M  Inca 

•  Moored  M  In  ea 

•  Close  Tethered  Mines 

•  Bottom  M  Ines 

•  B  u  rled  M  ines 

•  Rocket  Propelled  Mines 

•  In  ert  D  e  coy  Mines 

-  Non  -M  Ine  Bottom  Objects 


LAND  M  INES  : 

•  Anti-Personnel  Mines 

•  Anti-Tank  (Vehicle)  Mines 

•  Anti-Invasion  M  ines 


OBSTACLES 


OBSTACLES: 

-  Concertina  Wire 
-Anti-Tank  Tetrahedral 

-  Antl-T a  n  k  M  oats 

-  T  rip  W  ire  s 

-  Hedgehogs 
-Engineering  Stakes 
-Concrete  Blocks 

-  W  ooden  P  oles 

-  Booby  Traps 


Figure  1.  Littoral  Mine  Threat.  From  Rhodes,  1998. 


Environmental  data  collection  in  a  potential  adversaries’  littoral  region  is  often 
hampered  by  inaccessibility.  This  lack  of  accurate  data  causes  a  certain  degree  of 
uncertainty  in  MCM  planning.  As  a  result,  several  numerical  models  for  predicting  mine 
impact  burial  (IB)  have  been  developed.  The  most  promising  IB  model  was  originally 
developed  by  Amone  and  Bowen  in  1980.  Later  improvements  by  Satkowiak  (1987), 
Hurst  (1992)  and  others  have  resulted  in  the  current  version,  IMPACT  25.  IMPACT  25 
creates  a  two-dimensional  time  history  of  a  cylindrical  mine  as  it  falls  through  air,  water 
and  sediment  phases.  IB  prediction  is  largely  calculated  from  the  marine  sediment 
characteristics  and  mine  impact  orientation  and  velocity. 

Several  studies  with  regard  to  IB  have  been  conducted  over  the  years.  Taber 
(1999),  Smith  (2000)  and  Lott  (1995)  conducted  experiments  to  verify  the  IMPACT  25 
sediment  phase  calculations.  However,  accurate  IB  prediction  also  requires  that  the 
mine’s  trajectory  through  the  air  and  water  phases  is  adequately  accounted  for.  IMPACT 
25  approximates  mine  shape  as  a  cylinder,  includes  a  torque  adjustment  for  cases  where 
COM  does  not  coincide  with  the  center  of  buoyancy  and  calculates  the  air  and  water 
phase  trajectories  entirely  from  momentum  equations.  As  a  result,  calculated  trajectories 
tend  to  be  arc  shaped  and  lack  movement  of  the  mine  about  its  own  axis.  In  order  to 
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examine  what  effect  a  varying  center  of  mass  has  on  a  mine's  water  phase  trajectory,  the 
first  in  a  series  of  Mine  Drop  Experiments  (MIDEX)  was  conducted  in  June  200 1. 
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n.  MINE  DROP  EXPERIMENT  (MIDEX) 


A.  PREPARATION 

MIDEX  basically  consisted  of  dropping  each  of  three  right  cylinders  into  the 
water  where  each  drop  was  recorded  underwater  from  two  viewpoints.  The  controlled 
parameters  for  each  drop  were:  center  of  mass  position  (COM),  initial  velocity  (Vinit), 
drop  angle  and  the  ratio  of  mine’s  length  to  diameter.  Figure  2  depicts  the  overall  setup. 


Figure  2.  Equipment  used.  A  denotes  drop  angle  device,  B  mine  injector,  C  infrared  light 
sensor,  D  output  to  universal  counter,  E  mine  shapes. 

1.  Mine  Shapes  and  Center  of  Mass  Positions 

Three  mine  shapes  were  used  for  the  experiment.  All  had  a  circular  diameter  of  4 
cm,  however  the  lengths  were  15,  12  and  9  cm  respectively.  The  bodies  were  constmcted 
of  rigid  plastic  with  aluminum- capped  ends.  Inside  each  was  a  threaded  bolt,  running 
lengthwise  across  the  mine,  and  an  internal  weight  (Fig  3).  The  internal  weight  was  used 
to  vary  the  mine’s  COM  and  could  be  adjusted  fore  or  aft. 


Figure  3.  Internal  Components  of  Mine  Shape 
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COM  positions  were  denoted  using  2,  1,  0,  -1,-2.  COM  0  cases  were  identical  to 
the  IMPACT  25  model’s  uniform  density  assumption  (COM  =  CB).  All  other  cases 
indicated  the  relative  position  of  the  COM  to  the  CB.  Figure  4  displays  the  various 
measurements  for  each  COM  position. 


-  M 


MODEL  s  1 

L.1S13E9cm  D.4crn  m-2  7cn 
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Figure  4.  Mine  Shape  Characteristics.  The  leftmost  column  indicates  COM  position  0, 
while  the  rightmost  indicates  COM  position  2. 

2.  Model  Mine  Scaling 


Our  goal  was  to  choose  a  scale  that  was  somewhat  representative  of  the  real  world 
ratio  of  water  depth  to  mine  length,  but  at  the  same  time  would  be  large  enough  to  film 
and  would  not  damage  the  pool's  bottom.  The  model  mines  were  based  on  the  realistic 
assumption  that  a  3  m  mine  is  laid  in  water  depths  of  45  m,  thus  producing  a  15:1  ratio. 
This  ratio  was  close  to  our  baseline  15  cm  mine  being  dropped  into  a  pool  of  2.4  m  depth. 
The  addition  of  a  12  and  9  cm  length  allowed  for  later  comparison  of  the  sensitivity  of 
water  phase  trajectory  to  the  ratio  of  mine  length  over  diameter. 

3.  Initial  Velocity 


Initial  velocity  was  calculated  by  using  the  voltage  return  of  an  infrared  photo 
detector  located  at  the  base  of  the  mine  injector.  The  infrared  sensor  produced  a  square 
wave  pulse  when  no  light  was  detected  due  to  blockage  caused  by  the  mine’s  passage. 
The  length  of  the  square  wave  pulse  was  converted  into  time  by  using  a  universal 
counter.  Dividing  the  mine’s  length  by  the  universal  counter’s  time  yielded  Vinit.  The 
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mines  were  dropped  from  several  positions  within  the  injector  mechanism  in  order  to 
produce  a  range  of  Vinit. 

The  method  used  to  determine  Vinit  required  that  the  infrared  light  sensor  be 
located  above  the  water's  surface.  This  distance  was  held  fixed  throughout  the  experiment 
at  10  cm. 

4.  Drop  Angle 

Drop  angle  was  controlled  using  the  drop  angle  device.  Five  screw  positions 
marked  the  15,  30,  45,  60,  and  75-degree  positions.  The  drop  angles  were  determined 
from  the  lay  of  the  pool  walkway,  which  was  assumed  to  be  parallel  to  the  water’s 
surface. 

A  range  of  drop  angles  was  chosen  to  represent  the  various  entry  angles  that  air 
and  surface  laid  mines  exhibit.  This  range  produced  velocities  whose  horizontal  and 
vertical  components  varied  in  magnitude.  This  allowed  for  comparison  of  mine  trajectory 
sensitivity  with  the  varying  velocity  components. 

5.  Coordinate  System 

Two  coordinate  systems  are  used  to  describe  the  mine  falling  through  the  water 
column.  The  first  is  the  earth’s  coordinate  system  that  follows  the  right  hand  rule  (Fig  5) 
and  the  second  is  the  mine’s  coordinate  system  that  remains  rigidly  connected  to  the  mine 
(Fig  6).  Chapter  IV  provides  a  more  detailed  explanation  of  the  coordinate  systems  used. 


Figure  5.  Earth  Coordinate  System. 
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Figure  6.  Mine’s  Coordinate  System.  \|/  is  the  attitude  or  angle  in  the  z- plane 

Two  10  cm  grids  were  affixed  to  each  pool  wall.  Each  grid  was  constructed  out  of 
fiberglass  and  was  used  to  record  the  mine's  position  in  the  x,  -z  and  -y,  -z  planes  (Fig  7). 


Figure  7.  Background  grid.  Affixed  to  x,  -z  and  -y,  -z  planes 

B.  METHODOLOGY 

For  each  mn  the  mines  were  set  to  a  COM  position.  For  positive  COM  cases,  the 
mines  were  placed  into  the  injector  so  that  the  COM  was  located  below  the  center  of 
buoyancy.  For  negative  cases,  the  COM  was  located  above  the  center  of  buoyancy  prior 
to  release.  A  series  of  drops  were  then  conducted  in  order  of  decreasing  mine  length  for 
each  angle,  starting  at  15°  (Table  1).  Each  video  camera  had  a  film  time  of  approximately 
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one  hour.  At  the  end  of  the  day,  the  tapes  were  replayed  in  order  to  determine  clarity  and 
optimum  camera  position. 


15° 

o 

O 

CO 

45° 

60° 

75° 

mssssissnsMii 

2 

13 

15 

15 

15 

12 

1 

9 

15 

15 

15 

9 

0 

12 

15 

14 

18 

6 

-1 

6 

6 

6 

-2 

2 

6 

6 

0 

0 

Table  1.  Number  oi 

r  Drops  Conducted  by  Drop  Ang] 

le  and  COM  Position 
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HI.  DATA  RETREIVAL  AND  ANALYSIS 


A.  DATA  RETREIVAL 

Upon  completion  of  the  drop  phase,  the  video  from  each  camera  was  converted  to 
digital  format.  The  digital  video  for  each  view  was  then  analyzed  frame  by  frame  (30  Hz) 
in  order  to  determine  the  mine's  position  in  the  x,  -z  and  -y,  -z  planes.  The  mine's  top  and 
bottom  positions  were  input  into  a  MATLAB  generated  grid,  similar  to  the  ones  within 
the  pool.  The  first  point  to  impact  the  water  was  always  plotted  first.  This  facilitated 
tracking  of  the  initial  entry  point  throughout  the  water  column.  The  cameras  were  not 
time  synced;  thus,  the  first  recorded  position  corresponded  to  when  the  full  length  of  the 
mine  was  in  view. 

B.  SOURCES  OF  ERROR 

There  were  several  sources  of  error  that  hindered  the  determination  of  the  mine's 
exact  position  within  the  water  column.  Locations  above  or  below  the  camera's  focal 
point  were  subjected  to  parallax  distortion.  Placing  the  cameras  as  far  back  as  possible, 
while  still  being  able  to  resolve  the  individual  grid  squares,  rninimized  this  error.  Second, 
the  background  grids  were  located  behind  the  mine's  trajectory  plane.  This  resulted  in  the 
mine  appearing  larger  than  normal.  This  error  was  minimized  by  not  allowing  the  plotted 
points  to  exceed  the  particular  mine's  length.  Third,  an  object  injected  into  the  water  will 
generate  an  air  cavity.  This  air  cavity  can  greatly  affect  the  initial  motion,  particularly  at 
very  high  speeds  (hydro  ballistics).  The  air  cavity  effect  was  deemed  to  be  minimal  due 
to  the  low  inject  velocities  used. 

C.  DATA  ANALYSIS 

The  2-D  data  provided  by  each  camera  was  first  used  to  produce  raw  2-D  plots  of 
the  mine's  trajectory.  Next,  2-D  data  from  both  cameras  was  then  fused  to  produce  a  3-D 
history.  This  3-D  history  was  then  made  non-dimensional  in  order  to  generalize  the 
results.  The  non-dimensional  data  was  used  to  generate  impact  scatter  plots  and  was  also 
used  in  multiple  linear  regression  calculations.  The  non-dimensional  conversions  used 
were: 
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t  = 


dt 


V  T 

— ;COM= 
■v/gL  D 


2M  _  (x,y,z) 
L  ’  L 


,  u,\vw  ,  cos^rQp  angie) 

VsL 


(1) 


where  L  =  mine  length,  D  =  mine  diameter,  g  =  gravitational  force,  (x,y,z)  are  the  mine's 
earth  coordinates  and  (u,v,w)  are  the  nine's  velocities  in  the  x,  y  and  z  directions.  Table  2 


provides  a  summary  of  the  commonly  used  non-dimensional  COM  terms. 


Mine  Lenath  ► 

15 

12 

9 

COM  Position 

2 

0.1939 

0.1594 

0.1198 

1 

0.0969 

0.0797 

0.0599 

0 

0 

0 

0 

-1 

-0.0969 

-0.0797 

-0.0599 

-2 

-0.1939 

-0.1594 

-0.1198 

Table  2.  Commonly  used  non-dimensional  COM  positions.  Non-dimensional  COM 


positions  determined  using 


2M 

L 


,  where  M  is  given  in  Figure  4. 


1.  2-D  Plots 

The  raw  2-D  data  was  plotted  in  the  x,  -z  and  -y,  -z  planes  for  two  purposes.  First, 
it  provided  a  check  on  the  methodology  used  for  position  recording.  Second,  the  2-D  -y,  - 
z  data  allowed  us  to  develop  generalized  trajectory  patterns  that  the  mine's  appeared  to 
follow.  Figure  8  depicts  an  example  of  a  2-D  plot.  Trajectory  patterns  are  presented  in 
Chapter  5  and  Appendix  A. 

2.  Impact  Attitude 

The  mine's  final  data  point  was  considered  to  be  the  impact  point.  Impact  points 
were  grouped  by  COM  position  and  mine  length  in  order  to  establish  the  relationship  of 
COM  position  to  impact  attitude.  Figure  9  provides  a  summary  of  all  cases.  Further 
histograms  are  provided  in  Appendix  B. 

3.  Impact  Point  Locations 

The  lateral  movement  a  mine  experiences  after  drop  is  of  extreme  importance  to 
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Figure  9.  Impact  Angles  for  all  cases. 


the  MCM  planner.  When  the  threat  of  mines  exists,  a  mine  danger  area  is  established. 
The  greater  the  uncertainty  in  the  mine's  location,  the  larger  the  danger  area.  Increases  in 
a  mine  danger  area  translate  into  increased  clearance  times.  The  impact  points  were 
analyzed  to  determine  the  lateral  movement  experienced  as  influenced  by  COM  position, 
and  drop  angle  (Fig  10).  Additional  impact  plots  can  be  found  in  Appendix  C. 
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Impact  Points  (All  Cases) 


X  (non-dimensional) 

Figure  10.  Impact  Points  for  all  cases  and  conditions.  The  yellow  areas  represent  the 

cameras  used. 
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IV.  DYNAMICS  OF  A  CYLYNDRICAL  MINE 


A.  TWO  COORDINATE  SYSTEMS 

Two  coordinate  systems  are  used  to  describe  a  cylindrical  mine  falling  through 

the  water  column:  earth  and  body  (mine)  coordinates.  The  earth  coordinate  system  is 

fixed  to  the  swimming  pool  with  horizontal  sides  as  x  and  yaxis,  and  vertical  direction  as 
the  z-axis  (Fig.  5).  The  body  coordinate  is  rigidly  connected  with  the  cylindrical  mine. 

The  origin  of  the  body  coordinate  system  coincides  with  the  center  of  gravity  (CG);  the 

axis-ri  is  along  the  centerline  of  the  cylinder;  the  axis-r2  is  perpendicular  to  the  plane 
constructed  by  axes-ri  and  axis-z  {r\-z  plane);  and  the  axis-r3  lies  in  the  (n -z)  plane  and  is 
perpendicular  to  axis-ri.  The  selection  of  axes  (v,  y,  z.)  and  fa,  r2,  r3)  follows  the  right- 
hand  mle.  The  angles  that  the  three  axes  r\,  n,  r3  form  with  the  vertical  (upward  positive) 
are  called  (f>i,  ^2,  and  <f>3  (Fig.  11).  The  angle  c|)  1  is  usually  called  the  attitude  of  the  mine. 


Figure  11.  Mine  Coordinate  System.  The  axis  r3  originates  at  CG  and  comes  out  of  the 
page  towards  the  viewer.  CB  is  the  center  of  buoyancy  (volumetric  center)  of  the  mine. 

B.  HYDRODYNAMIC  THEORY  OF  MINE  IMPACT  BURIAL 

The  essential  elements  of  the  mine  impact  burial  model  translate  into  the  science 
and  engineering  of  hydrodynamic  processes  associated  with  a  falling  object  and  sediment 
transport.  Any  solid  object  fading  through  a  fluid  (air  and  water)  should  obey  two 
physical  principles:  (2)  momentum  balance  and  (3)  moment  balance: 


15 


n(dV*/dt*)  dm*  =  W*+Fb*+  Fd* 

(2) 

n[r*x(dV*/dt*)J  dm*  =M*, 

(3) 

where  the  superscript  *  denotes  dimensional  variables.  V  is  the  velocity  of  the  mine,  W 
the  gravitational  force,  Fb  the  buoyancy  force,  Fd  the  drag  force,  and  M  the  resultant 
moment  of  momentum.  Existing  IBPM  models  only  consider  the  momentum  balance  of 
the  mine  [i.e.,  Eq.  (2)]  and  disregard  the  moment  balance  of  the  mine  [i.e.,  Eq.  (3)].  Such 
an  incomplete  hydrodynamics  of  the  model  leads  to  unrealistic  prediction  of  the  mine's 
water  phase  trajectory  (Fig.  12).  By  considering  momentum  and  moment  balance,  the 
falling  object  should  have  a  spiral-type  motion  (Fig  13).  Without  the  spiral-type  motion, 
the  IBPM  may  over-predict  the  impact  burial  depth  (Chu  et  al.  2000). 
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Figure  12.  Mine  Motion  without 


With  Constant  Rotnbon 


Consideration  of  Moment  Equation. 
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With  Moment  Equations 
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Spiral  -Type  Motion  ^ 

|  (I 

Figure  13.  Mine  Motion  with  Consideration  of  both  Momentum  and  Moment  Equations. 


Let  (V*X,V*)  be  the  three  components  of  the  velocity  of  CG  and  ((Oj  ,or2.OJ. ) 
the  components  of  the  angular  velocity,  referring  to  the  direction  of  the  mine- fixed 
coordinate  system.  The  momentum  equation  (2)  becomes: 


dV T 
dt 


+co2y3  -co3y. 


: <?p-  p-  cos¥,  i  (v;  -yj,) 

P,„  2Lpffl 


(4a) 


+w3v;  -co,v;  =  g  Pjh — p«i  cosy ,  \v*  i  cy2*  -y;2)  (4b) 

dt  pm  “  2Lp,„ 


dVs-  +co y;  -co2v; = g  Pm  Pw  cost)/ 3  - i  v  \  (y;  - y;3 )  (4c) 


dt 


2^P„ 


where  Cd  is  the  drag  coefficient,  g  the  gravitational  acceleration,  L  the  length  of  the  mine, 
and  (Kw  ,yj2>yj3)  is  the  water  velocity.  The  independent  and  dependent  variables  are 
made  non-dimensional  by: 


t 


The  momentum  equations  (4a)- (4c)  become: 
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dV, 


C  0  IV  0  —  D 

_  D  rw  I  1/  I  f\7  T /  \  i  r  m  rM 


-^+0)^3 -co3V2  =  -^1 V I  (M -V*)  + 

*  2p„,  p 


C0SV| /, 


Co  w  0  —  D 

_  11/1/17  1/  \  I  r  m  rv 


^+mjV,-0),V3  =  -^IVI(V!-V„j)  + 

*  2pm  P 


COSt|/, 


^+0),v2  — co2V,  =  V  |  (v3  -  V  )  +  Pm  Ph  cosy, 

*  -  "  2p„,  '  pm 


The  non-dimensional  equation  of  the  moment  of  momentum  (3)  becomes: 


j  d C0i  ,  j  ,  ,  .dco3  ,  LM* 

/, -  +  {J3  -  J2) CO,C03  -  /31( i+CO^)  = L 

dt  dt  ”  g 


J 2  mi+  (7(  _  /3)C03C01  -  /31  (C03  —  COf )  = 

dt  ~~  g 


T  d CO,  .  7  _ .  ,  .dco,  ,  LM* 

7 3  — 1 +  (J2-J1 )co1co2  — /31( — 1 — co2co3 )  = - -, 

dt  dt  g 

where  J\ ,  J2,  and  J3  are  the  three  moments  of  gyration, 

7,  =| (r2  +  r2)dm*  ,  J2  =  J (r2  +  r2)dm  ,  J3  =  J (r2  +  r2)dm 

and  the  moment  of  deviation  (or  inertia  products  of  second  order), 

J,  |  =  |  rfflm 

The  orientation  of  the  mine  (()  1 ,  p2-  ps)  is  determined  by 


d 

— COSl|/,  =co3  COS\|/2 -C02  COS\|/3 

dt 


d 

— cost)/ 2  =00,  COSt|/3  —  C03  cost)/! 

dt 


d 

— COSt)/3  =  C0^  COSt|/,  -co,  cost)/2 . 

dt 


The  nine  non-dimensional  equations  (6a-c),  (7a-c),  and  (lOa-c)  are  the 
describing  the  mine  movement  in  the  water  column. 


(6a) 

(6b) 

(6c) 

(7a) 

(7b) 

(7c) 

(8) 

(9) 

(10a) 

(10b) 

(10c) 

basic  system  for 
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C.  PHYSICAL  CHARACTERISTICS  OF  MINE 

For  a  cylindrical  mine,  two  characteristics  are  important  for  the  prediction  of  mine 
movement  in  the  water  column:  the  ratio  between  length  and  diameter  (L/D)  and  the  ratio 
M/2L  with  M  being  the  distance  between  the  centers  of  gravity  and  buoyancy.  Positive 
(negative)  values  of  M  refers  to  the  center  of  gravity  below  (above)  the  center  of 
buoyancy.  In  MIDEX,  L/D  =  15/4,  12/4,  and  9/4;  M/2L  values  are  provided  in  Table  2. 

D.  INITIAL  CONDITIONS  IN  MIDEX 

During  the  experiment,  the  initial  conditions  were  initial  velocity  and  entry  angle. 
Thus,  for  the  dynamical  system  (6a-c),  (7a-c),  and  (lOa-c),  the  initial  conditions  at  t  =  0 
were: 

Vj=-^=,  V2  =  0,  V3  =  0,  Yi  =  y2  =  y3  =  0,  ^  =  Drop  Angle,  <|)2  =  <b  =  0.  (11) 
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V.  RESULTS 


A.  TRAJECTORY  PATTERNS 

By  analyzing  the  2-D  -y,  -z  planar  plots  we  were  able  to  develop  seven  general 
trajectory  patterns.  Pattern  names  are  based  upon  the  mine  coordinate  system.  The  -y,  -z 
planar  plots  were  chosen  for  trajectory  analysis,  as  this  plane  was  parallel  to  the  direction 


of  mine  drop.  The  generalized  trajectory  patterns  are  described  in  Table  3  and  Figures 
14a  and  b.  Appendix  A  contains  all  of  the  remaining  analyzed  2-D  plots. 


Mine  Trajectory  Pattern 

Description 

Straight  or  Slant 

Mine  exhibited  little  angular  change  about  z-axis.  For  straight 
mine  attitude  remained  nearly  parallel  with  z-axis  (±  15°).  For 
slant,  mine  attitude  was  45°  off  z-axis  (±  15°). 

Spiral 

Mine  experienced  rotation  about  z-axis  throughout  its  water 
phase  trajectory. 

Flip 

Initial  water  entry  point  rotated  at  least  180°  during  mine 
motion. 

Flat 

Mine's  angle  with  vertical  near  90°  for  most  of  the  trajectory. 

Seesaw 

Similar  to  the  flat  pattern  except  that  mine's  angle  with  vertical 
would  oscillate  between  greater  (less)  than  90°  and  less 
(greater)  than  90°  -  like  a  seesaw. 

Combination 

Complex  trajectory  where  mine  exhibited  several  of  the  above 
patterns. 

Table  3.  Description  of  Mine  Coordinate  Based  Trajectory  Patterns 

B.  IMPACT  ATTITUDE 


IB  is  largely  determined  from  the  impact  attitude  of  the  mine.  Mines  whose 
impact  attitudes  are  perpendicular  (vFi  =  0  or  180°)  to  the  sediment  interface  will 
experience  the  largest  degree  of  IB.  (Taber  1999).  It  is  therefore  important  to  analyze  the 
relationship  between  impact  attitude  and  the  controlled  parameters,  drop  angle,  Ymt-  L/D 
and  COM  position.  Both  L/D  and  Vmit  had  little  influence  on  impact  attitude.  COM 
position  and  drop  angle,  however,  were  the  largest  determinants  of  mine  impact  attitude. 

From  Figure  9  it  is  apparent  that  there  are  several  peaks  centered  near  90°,  140° 
and  180°.  Further  analysis  reveals  that  these  peaks  correspond  to  COM  positions  0,  1  and 
2  respectively  (Fig.  15).  COM  positions  -1  and  -2  followed  the  same  trend  as  their 
positive  counterparts.  In  our  coordinate  system  the  attitude  is  measured  with  respect  to 
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the  mine's  first  entry  point.  When  looking  into  the  -y,  -z  plane,  a  mine  with  a  CP  or  40° 
attitude  will  look  the  same  as  a  mine  with  a  180°  or  140°  attitude. 

Although  drop  angle  was  not  the  most  influential  parameter,  variations  did  induce 
changes  in  impact  orientation.  As  drop  angle  increased,  the  likelihood  of  any  lateral 
movement  decreased.  This  allowed  for  impact  angles  that  were  more  vertically 
orientated.  This  is  primarily  due  to  the  fact  that  the  vertical  components  of  velocity  were 
greater  than  those  at  shallow  angles.  Thus,  the  time  to  bottom  and  time  for  trajectory 
alteration  was  less. 


Straight  Spiral 

L=  15cm;  Vi=  1  835m/S  L=  9cm;  Vi=  1  553m/s 
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E  -1 


-1.5 


Flat 

L=  15cm;  Vi=  2  18m/s 


-05 


Y(m> 


Y(m) 


Y  (m) 


Figure  14a.  Trajectory  Examples. 
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Flip 


Seesaw 


Flip -Flat 


L=15cm:  Vi=  2  874m/s 


L=  15cm:  Vi=  1  379m/s 


L=  9cm:  Vi=  2  294m/s 


Y<m>  Y(m)  Y(m) 


Figure  14b.  Trajectory  Examples 


Impact  /-not  Freqwney  of  Ocsjrranee 


Itcbcj  An^c  Oroiccd  in  5*  0m 

Figure  15.  Relationship  between  COM  Position  and  Impact  Attitude. 
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c. 


IMPACT  POINTS 


Using  a  methodology  similar  to  that  used  in  impact  attitude  analysis,  impact  point 
scatter  plots  were  analyzed  by  the  controlled  parameters,  drop  angle,  Vinit  and  COM 
position.  Impact  point  scatterplots  are  provided  in  Appendix  C.  COM  2  and  0  cases  fell 
near  the  drop  point  greater  than  90%  of  the  time  while  COM  1  cases  displayed  the  most 
variability.  Additionally,  the  flip  experienced  in  negative  COM  cases  induces  a  greater 
degree  of  lateral  movement  than  in  positive  COM  cases. 

D.  MULTIPLE  LINEAR  REGRESSION  ANALYSIS 

Multiple  linear  regression  analysis  was  performed  on  the  last  recorded  points 
(impact  points)  using  a  least  squares  technique  as  outlined  by  Walpole  (1998).  The 
purpose  of  this  was  to  establish  a  relationship  between  the  input  non-dimensional 
parameters;  drop  angle,  COM  position,  Vmit  and  L/D,  and  output  variables;  (Xm,ym), 
(u,v,w)  and  V|/ 1  (Table  4).  The  basic  equation  used  was: 

f  =B0  +6^  +B2x^  +B3x3  +B4x4  +  e;  (12) 


where  f,  is  the  desired  output  parameter,  (Xm,  ym,  i|/|.  u,  v,  w),  betas  are  the  correlation 
coefficients,  xi  is  drop  angle,  x2  is  L/D,  x3  is  Vinit  and  X4  is  COM  position.  The  results 


Xm 

ym 

¥1 

u 

V 

w 

Po 

-.0746 

-.0546 

102.5691 

.0040 

-.0135 

-.9481 

pi 

.1190 

-.0828 

-13.3508 

-.0075 

-.0106 

-.1080 

P2 

-.0469 

-.0798 

-.5009 

-.0011 

.0005 

.0295 

P3 

.0372 

.0622 

1.0437 

.0025 

.0011 

-.0221 

P4 

.2369 

.4330 

472.2135 

-.0090 

.0537 

-1.2467 

Table  4.  Multiple  Linear  Regression  Analysis  Correlation  Coefficients. 


indicate  that  COM  position  has  the  largest  influence  on  all  output  variables.  As  a  check 
we  determined  \|/i  for  a  case  where  L/D  =  3.75,  Drop  angle  =  15°,  COM  =  .1939  and  Vinit 

=  3m/s.  After  non-dimensionalizing  drop  angle  and  Vinit  we  found  that  the  impact  attitude 
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was  181.2°,  in  good  agreement  with  observations.  Similarly,  for  COM  =  .0969  and  0  we 
found  that  the  impact  attitudes  were  136.1°  and  90.4°  respectively.  This  is  also  in  good 
agreement  with  observations  and  follows  the  trends  established  in  the  histogram  analysis. 
Impact  point  data  tables  can  be  found  in  Appendix  D. 
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VL  DISCUSSION 


The  water  phase  trajectory  a  mine  experiences  ultimately  determines  the  impact 
orientation.  In  MIDEX,  the  categorizing  of  trajectories  into  general  patterns  served  two 
purposes.  First,  the  pattern  name  gives  a  sense  of  impact  attitude.  Second,  it  allowed  for  a 
visualization  of  the  motion  a  mine  experiences  about  its  own  axis.  Observed  trajectories 
were  found  to  be  most  sensitive  to  COM  position,  drop  angle  and  mine  length.  As  COM 
distance  increased  from  the  CB  the  mine  tended  to  follow  a  straight  pattern.  As  COM  was 
moved  closer  to  the  CB  the  mine's  trajectory  tended  towards  being  more  parallel  with  the 
pool's  bottom  (Tables  5a-d).  At  steep  drop  angles,  the  mine  experienced  little  lateral 
movement  and  tended  towards  a  straight  pattern.  Additionally,  as  mine  length  decreased 
more  complex  trajectory  patterns  developed.  This  included  significant  rotation  about  the 
vertical  axis  and  increased  lateral  movement. 


15 

12 

9 

DroD  Anale:  15° 

Straiaht 

Straiaht 

Straiaht-Slant 

Slant-Straiaht 

SDiral 

SDiral 

DroD  Anale:  30° 

Slant-Straiaht 

Slant 

SDiral 

Slant-Straiaht 

Straiaht 

SDiral 

Slant-Straiaht 

Slant-Straiaht 

SDiral 

Straiaht-SDiral 

Straiaht-SDiral 

SDiral 

Straiaht 

Straiaht 

Spiral 

Slant 

Slant-Straiaht 

Slant-SDiral 

DroD  Anale:  60° 

Straiaht 

Straiaht-Slant 

SDiral 

Straiaht 

Straiaht 

SDiral 

Droo  Anale:  75  ° 

Straiaht 

Straiaht 

Spiral 

Straiaht 

Straiaht-SDiral 

Slant 

Straiaht 

Straiaht 

SDiral 

_ Straight _ 

_ Straiaht-Spiral _ 

_ Straiaht-Spiral _ 

Table  5a.  Observed  Trajectory  Patterns  for  COM  Position  2. 
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COM  Position:  1 

Mine  Lenath: 

15 

12 

9 

DroD  Anale:  15° 

Slant 

SDiral 

Straiaht-SDiral 

DroD  Anale:  30° 

Slant 

Slant-SDiral 

SDiral 

Straiaht-Slant 

Slant 

SDiral 

DroD  Anale:  45° 

Slant 

Slant 

SDiral 

Straiaht 

Straiaht-SDiral 

Slant 

Slant 

Straiaht 

Straiaht-SDiral 

Slant 

Straiaht 

Straiaht-SDiral 

DroD  Anale:  75  ° 

Straiaht 

Straiaht 

Straiaht-SDiral 

Straiaht 

Straiaht 

Straiaht-SDiral 

Straight 

Straight-Spiral 

Slant-Spiral 

Table  5b.  Observed  Trajectory  Patterns  for  COM  Position 


15 

12 

9 

DroD  Anale:  15° 

Seesaw 

Flat 

Flat-Seesaw 

Flat 

DroD  Anale:  30° 

Seesaw 

Flat 

SDiral-Seesaw 

Flat 

Flat 

Flat 

Flat-Soiral 

Straiaht-Flat 

Straiaht-Flat 

Straiaht 

Straiaht-Flat 

Straiaht-Flat-SDiral 

DroD  Anale:  60° 

SDiral-Seesaw 

Straiaht-Seesaw 

Straiaht-Flat 

Straiaht-Flat-SDiral 

Straiaht 

Straiaht-Seesaw-SDiral 

Straiaht-SDiral-Seesaw 

Straiaht-SDiral-Seesaw 

Straiaht-Seesaw 

Straight-Spiral-Seesaw 

_ Straight _ 

_ Straight-Flat-Spiral _ 

Table  5c.  Observed  Trajectory  Patterns  for  COM  Position  0. 


15 

12 

9 

DroD  Anale:  30° 

FliD-Straiaht 

FliD-Slant 

FliD-Straiaht-SDiral 

DroD  Anale:  45° 

FliD-Straiaht 

FliD-Straiaht-SDiral 

FliD-Straiaht 

FliD-Straiaht 

Flip-Straight 

FliD-Straiaht-SDiral 

DroD  Anale:  30° 

FliD-Straiaht 

FliD-Slant 

Straiaht-FliD-Seesaw 

DroD  Anale:  60° 

FliD-Straiaht 

FliD-Straiaht 

FliD-SDiral-Seesaw 

Straight-Flip 

Slant-Flip-Slant 

Flip-Spiral 

Tab: 


e  5d.  Observed  Trajectory  Patterns  for  Negative  COM  Cases. 


In  order  to  evaluate  the  validity  of  the  IMPACT  25  water  phase  trajectory 
patterns,  a  comparison  between  our  results  and  the  model's  was  performed.  The  mine 
shape  characteristics  (length,  diameter,  COM  distance),  drop  angle,  release  altitude  (.1 
m),  water  temperature  (23.8°C)  and  initial  velocities  (converted  into  horizontal  and 
vertical  components)  were  entered  into  IMPACT  25.  Model  output  was  organized  by 

COM  position,  mine  length,  drop  angle  and  impact  angle  (Table  6). 
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coordinate  system  used  during  MIDFX. 


Overall,  the  IMPACT  25  model  decreased  the  impact  angle  as  COM  distance 
decreased.  This  agreed  with  the  observed  trend  in  the  histogram  analysis.  Additionally, 
the  model  predicted  impact  angle  for  COM  2  cases  and  the  flip  for  negative  cases  fairly 
well.  However,  significant  deviation  between  observed  and  model  output  occurred  for 
COM  cases  1  and  0.  This  deviation  was  primarily  due  to  the  fact  that  the  trajectories  are 
calculated  using  only  momentum  equations.  This  caused  the  model  to  be  more  sensitive 
to  drop  angle,  which  resulted  in  the  impact  angles  being  more  vertical  as  drop  angle 
increased.  As  such,  the  model  tended  to  be  slow  in  predicting  alterations  in  trajectory 
caused  by  moving  the  COM  closer  to  the  CB.  Furthermore,  the  IMPACT  25  patterns 
tended  to  be  more  arc  shaped  and  did  not  include  the  spiraling  motion  frequently 
observed  for  the  9  cm  mine. 
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vn.  CONCLUSIONS 


MTDKX  is  the  first  step  in  an  ongoing  process  to  better  understand  and  predict  the 
various  parameters  that  effect  a  mine's  motion.  Although  the  cmdeness  of  the  experiment 
precluded  development  of  prediction  equations,  it  did  prove  one  key  point.  COM  position 
has  the  largest  influence  on  mine  water- phase  trajectories. 

The  observed  trajectories  were  far  more  complex  than  those  theorized  by  using 
only  the  momentum  and  rotation  equations.  Simply  entering  the  initial  attitude  and 
rotation  rate  into  the  model  will  not  satisfy  the  needs  of  the  modem  naval  warfighter.  At  a 
minimum,  updates  to  the  IMPACT  25  model  should  include  the  more  realistic  moment 
equations. 

Further  research  on  mine  hydrodynamics  is  needed.  The  research  needs  to  expand 
beyond  the  simple  cylindrical  shaped  mine  to  those  that  are  irregularly  shaped  (Rockan 
and  Manta  types).  Additionally,  the  utilization  of  scaled  down  versions  should  be 
explored.  A  smaller  mine  that  can  be  modeled  as  accurately  as  its  real  counterpart  will 
save  time,  money  and  will  require  less  logistical  support. 
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APPENDIX  A.  2-D  PLOTS 


Appendix  A  contains  all  of  the  2-D  -y,  -z  plots  that  were  used  to  establish 
generalized  trajectory  patterns.  The  scales  have  been  left  in  dimensional  units  and  the 
camera  viewpoint  is  perpendicular  to  the  mine's  drop  direction  (-y  axis).  The  pattern 
names  were  chosen  based  upon  the  general  shape  of  the  trajectory  pattern  in  mine 
coordinates  (straight,  slant,  spiral,  flip,  flat,  seesaw  or  combination). 


L=  15cm;  Vi=  2.302m/s 


Drop  Angle  15;  COM:2 
L=  12cm;  Vi=  2  761  m/s 


L=  9cm;  Vi=  2.685m/s 
Straight-Slant 

•  / 


_l - 1 - 1 - 1 - 1_ 


Y(m) 


Y<m) 


Y(m) 
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Z(m) 


L=  15cm;  Vi=1.186m/s 


Drop  Angle:  15;  COM:2 
L=  12cm;  Vi=  1.406m/s 


L=  Scm;  Vi=  1.512m/s 


Y  (m) 


Y  (m) 


Y  (m) 
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Z(m) 


Drop  Angle:  15;  COM  2 

L=  15cm;  Vi=  2.1 9m/s  L=  12cm;  Vi=  1.91 8m/s  L=  9cm;  Vi=  2.397m/s 


Y  (m)  Y(m)  Y(m) 
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Z(m) 


Drop  Angle:  15;  COM  2 


L=  1 5cm;  Vi=  1 ,897m/s  L=  1 2cm;  Vi=  1 ,448m/s 


-0.6-0. 4-0.2  0  0.2  0.4  -0.5  0  0.5 

Y  (m)  Y(m) 


L=  9cm;  Vi=  1 ,72m/s 


Y  (m) 
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(w)  Z 


Drop  Angle:  30;  COM:2 

L=  15cm;  Vi=  1.1 53m/s  L=  12cm;  Vi=  1 ,259m/s  L=  9cm;  Vi=  1  122m/s 


Y  (m)  Y  (m)  Y  (m) 


37 


(w)  Z 


Drop  Angle:  30;  COM:2 

L=  1 5cm;  Vi=  2.082m/s  L=  1 2cm;  Vi=  1 ,953m/s  L 


-0.6-0. 4-0. 2  0  0.2  0.4  -0.6-0.4-0.2  0  0.2  0.4  -0.5 

Y  (m)  Y  (m) 


9cm;  Vi=  2.199m/s 


_ I _ I 

0  0.5 

Y(m) 
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(w)  Z 


L=  15cm;  Vi=  1.525m/s 


Drop  Angle:  30;  COM:2 
L=  12cm;  Vi=  1.791  m/s 


L=  9cm;  Vi=  1 ,562m/s 


Y(m) 


-0.5  0  0.5  -0.5  0 

Y  (m)  Y  (m) 
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(w)  Z 


L=  15cm;  Vi=  1.244m/s 


Y  (m) 


Drop  Angle:  30;  COM:2 

L=  9cm;  Vi=  1.26m/s 


-0.5  0  0.5  -0.5  0 

Y  (m)  Y  (m) 


L=  12cm;  Vi=  1.527m/s 


0.5 


40 


(w)  Z 


L=  15cm;  Vi=  2.309m/s 


Y  (m) 


Drop  Angle:  30;  COM:2 


L=  12cm;  Vi=2.724m/s 


L=  9cm;  Vi=  2.746m/s 


Y(m) 
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(w)  Z 


Drop  Angle:  45;  COM:2 
L=  1 5cm;  Vi=  1 ,328m/s  L=  1 2cm;  Vi=  1 ,565m/s 


L=  9cm;  Vi=  1 ,626m/s 


42 


(w)  Z 


Drop  Angle:  45;  COM:2 


L=  1 5cm;  Vi=  2.684m/s  L=  1 2cm;  Vi=  2.649m/s  L=  9cm;  Vi=  2.508m/s 


Y  (m)  Y  (m)  Y  (m) 
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(w)  Z 


Drop  Angle:  45;  COM:2 

L=  15cm;  Vi=  2.23m/s  L=  12cm;  Vi=  2.3m/s 


Y  (m)  Y  (m) 


L=  9cm;  Vi=  2.242m/s 

- 1 - 1 - 1 - 1 - 1 - 


_ I _ I _ I _ I _ I _ 

-0.6-0. 4-0. 2  0  0.2  0.4 
Y(m) 


44 


(w)  Z 


Drop  Angle:  45;  COM:2 


L=  1 5cm;  Vi=  1 ,832m/s  L=  1 2cm;  Vi=  1 ,836m/s  L=  9cm;  Vi=  1 ,823m/s 


45 


(w)  Z 


Drop  Angle:  45;  COM:2 

L=  15cm;  Vi=  2.891  m/s  L=  12cm;  Vi=  3.1 23m/s  L=  9cm;  Vi=  2.953m/s 


-0.5  0  0.5  -0.5  0  0.5  -0.8-0.6-0.4-0.2  0  0.2 

Y  (m)  Y  (m)  Y  (m) 
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(w)  Z 


Drop  Angle:  60;  COM:2 

L=  1 5cm;  Vi=  1 ,567m/s  L=  1 2cm;  Vi=  1 ,708m/s  L=  9cm;  Vi=  1 ,838m/s 
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(w)  Z 


Drop  Angle:  60;  COM:2 


L=  1 5cm;  Vi=  2.935m/s  L=  1 2cm;  Vi=  2.963m/s 


L=  9cm;  Vi=  3.221  m/s 


Y(m) 


48 


(w)  Z 


Drop  Angle:  60;  COM:2 

L=  1 5cm;  Vi=  2.709m/s  L=  1 2cm;  Vi=  2.546m/s  L= 


9cm;  Vi=  2.51 2m/s 

- 1 - r 


/ 

_ i _ i 

0  0.5 

Y(m) 


49 


(w)  Z 


Drop  Angle:  60;  COM:2 


L=  1 5cm;  Vi=  1 ,835m/s  L=  1 2cm;  Vi=  1 ,722m/s 


L=  9cm;  Vi=  1 ,942m/s 


Y(m) 


50 


(w)  Z 


Drop  Angle:  60;  COM:2 


L=  15cm;  Vi=  3.4m/s  L=  12cm;  Vi=  3.544m/s 


L=  9cm;  Vi=  3.607m/s 

- 1 - 1 - 


_ I _ I _ U 

-0.5  0  0.5 

Y(m) 
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(w)  Z 


Drop  Angle:  75;  COM:2 


L=  1 5cm;  Vi=  2.41 3m/s  L=  1 2cm;  Vi=  2.31 8m/s 


L=  9cm;  Vi=  2.25m/s 


-0.8-0. 6-0. 4-0. 2  0  0.2 


Y(m) 


52 


(w)  Z 


L=  15cm;  Vi=  3.31 5m/s 


Drop  Angle:  75;  COM:2 


L=  1 2cm;  Vi=  3.209m/s  L=  9cm;  Vi=  3.244m/s 


Y  (m)  Y  (m) 


53 


(w)  Z 


Drop  Angle:  75;  COM:2 


L=  1 5cm;  Vi=  2.671  m/s  L=  1 2cm;  Vi=  2.41 7m/s 


L=  9cm;  Vi=  2.496m/s 


Y(m) 
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(w)  Z 


Drop  Angle:  75;  COM:2 


L=  1 5cm;  Vi=  3.647m/s  L=  1 2cm;  Vi=  3.761  m/s 


L=  9cm;  Vi=  3.756m/s 


Y(m) 
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(w)  Z 


Drop  Angle:  15;  COM:1 


L=  15cm;  Vi=  1.04m/s  L=  12cm;  Vi=  1.11 8m/s 


Y  (m)  Y  (m) 


L=  9cm;  Vi=  1.51  m/s 


i  i  r 

Straight-Spiral 

- T 

■ 

' 

_ 1 _ 

1 

- 

_ 1 _ L 

-0.5  0  0.5 

Y(m) 
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(w)  Z 


Drop  Angle:  15;  COM:1 


L=  1 5cm;  Vi=  1 ,059m/s  L=  1 2cm;  Vi=  1 .21 2m/s  L=  9cm;  Vi=  1 ,709m/s 


Y  (m)  Y  (m)  Y  (m) 
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(w)  Z 


Drop  Angle:  15;  COM:1 


L=  15cm;  Vi=  2.332m/s 


Y  (m) 


L=  12cm;  Vi=  1.894m/s  l 


Y(m) 


9cm;  Vi=  1.88m/s 


0  0.5 

Y(m) 
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(w)  Z 


Drop  Angle:  30;  COM:1 


L=  1 5cm;  Vi=  0.91 7m/s  L=  1 2cm;  Vi=  1 ,238m/s 


-0.8-0. 6-0. 4-0. 2  0  0.2-0. 8-0. 6-0. 4-0. 2  0  0.2 
Y  (m)  Y  (m) 


L=  9cm;  Vi=  1 ,204m/s 


Y(m) 
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(w)  Z 


L=  15cm;  Vi=  1 ,92m/s 


Y  (m) 


Drop  Angle:  30;  COM:1 


L=  12cm;  Vi=  1 ,795m/s  L=  9cm;  Vi=  1 ,984m/s 


Y  (m)  Y  (m) 


60 


(w)  Z 


Drop  Angle:  30;  COM:1 


L=  1 5cm;  Vi=  1 .871  m/s  L=  1 2cm;  Vi=  1 ,778m/s 


Y  (m)  Y  (m) 


L=  9cm;  Vi=  1.771  m/s 


-0.6-0. 4-0. 2  0  0.2  0.4 
Y(m) 
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(w)  Z 


Drop  Angle:  30;  COM:1 
L=  12cm;  Vi=  1.604m/s 


L=  15cm;  Vi=  1.347m/s 


L=  9cm;  Vi=  1 ,586m/s 


Y(m) 


Y(m) 


Y(m) 
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(w)  Z 


L=  15cm;  Vi=  3.079m/s 


Drop  Angle:  30;  COM:1 
L=  12cm;  Vi=  2.71 9m/s  L= 


9cm;  Vi=  2.98m/s 


i - r 


Y(m) 


63 


(w)  Z 


Drop  Angle:  45;  COM:1 
L=  12cm;  Vi=  1.483m/s 


L=  15cm;  Vi=  1 ,46m/s 


L=  9cm;  Vi=  1 ,553m/s 


Y  (m)  Y  (m)  Y  (m) 
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(w)  Z 


Drop  Angle:  45;  COM:1 

L=  1 5cm;  Vi=  2.623m/s  L=  1 2cm;  Vi=  2.788m/s  L=  9cm;  Vi=  2.74m/s 


Y  (m)  Y  (m)  Y  (m) 
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(w)  Z 


L=  15cm;  Vi=  2.203m/s 


Y  (m) 


Drop  Angle:  45;  COM:1 

L=  1 2cm;  Vi=  2.305m/s  L=  9cm;  Vi=  2.397m/s 


Y(m) 


Y(m) 


66 


(w)  Z 


Drop  Angle:  45;  COM:1 
L=  1 5cm;  Vi=  1 ,703m/s  L=  1 2cm;  Vi=  2.093m/s 


L=  9cm;  Vi=  2.033m/s 


Y  (m)  Y  (m) 


Y(m) 
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(w)  Z 


Drop  Angle:  45;  COM:1 


L=  1 5cm;  Vi=  2.945m/s  L=  1 2cm;  Vi=  2.803m/s 


L=  9cm;  Vi=  3.128m/s 


-0.6-0. 4-0. 2  0  0.2  0.4 


Y(m) 
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(w)  Z 


Drop  Angle:  60;  COM:1 


L=  1 5cm;  Vi=  1 ,8m/s  L=  1 2cm;  Vi=  1 ,865m/s 


Y  (m)  Y  (m) 


L=  9cm;  Vi=  1 ,946m/s 


Y(m) 
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(w)  Z 


L=  15cm;  Vi=  2.984m/s 


Drop  Angle:  60;  COM:1 

L=  1 2cm;  Vi=  3.497m/s  L=  9cm;  Vi=  3.289m/s 


Y(m) 


Y(m) 


70 


(w)  Z 


L=  15cm;  Vi=  2.91 2m/s 


Drop  Angle:  60;  COM:1 
L=  12cm;  Vi=2.929m/s 


L=  9cm;  Vi=  2.932m/s 


Y(m) 


Y(m) 


71 


(w)  Z 


Drop  Angle:  60;  COM:1 

L=  1 5cm;  Vi=  2.41 5m/s  L=  1 2cm;  Vi=  2.494m/s  L=  9cm;  Vi=  2.479m/s 


Straight-Spiral 


Y(m) 
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Drop  Angle:  75;  COM:1 

L=  1 5cm;  Vi=  2.358m/s  L=  1 2cm;  Vi=  2.484m/s  L=  9cm;  Vi=  2.428m/s 
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Drop  Angle:  30;  COM:0 
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Drop  Angle:  60;  COM:0 
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Drop  Angle:  45;  COM:-2 
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APPENDIX  B.  HISTOGRAM  PLOTS 


Appendix  B  contains  all  of  the  histogram  plots  that  were  derived  from  the  last 
recorded  data  point.  The  data  tables  used  for  histogram  production  are  found  in  Appendix 
D.  The  histograms  are  broken  down  by  mine  length  and  non-dimensional  COM  position. 
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impact  Ange  Frequency  for  COM  *  0 .159;  Length  =  12cm 
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Impact  Ange  Frequency  tor  COM  *  0  06;  Length  «  9cm 
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APPENDIX  C.  IMPACT  POINTS 


Appendix  C  contains  all  of  the  x-y  plane  impact  points.  The  impact  points  are 
derived  from  the  mine's  volumetric  center  x  and  y  coordinates  and  is  presented  in  three 
columns  and  five  rows.  The  columns  correspond  to  mine  length,  while  the  rows  represent 
the  COM  positions  beginning  with  position  2  to  position  -2.  There  are  two  sets  of  plots, 
one  for  all  drop  angles  and  another  set  for  each  drop  angle.  Each  of  the  individual  plots  is 
assigned  a  symbol  based  upon  the  mine's  Mnt-  The  velocity  ranges  used  are  in  m/s  and  are 
provided  in  the  legend  below: 
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Impact  Points  for  L/D  =  3.75  Impact  Points  for  L/D  =  3  Impact  Points  for  L/D  =  2.25 
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Impact  Points  for  L/D  =  3.75  Impact  Points  for  L/D  =  3  Impact  Points  for  L/D  =  2.25 
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Impact  Points  for  L/D  =  3.75  Impact  Points  for  L/D  =  3  Impact  Points  for  L/D  =  2.25 
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Impact  Points  for  L/D  =  3.75  Impact  Points  for  L/D  =  3  Impact  Points  for  L/D  =  2.25 
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Impact  Points  for  L/D  =  3.75  Impact  Points  for  L/D  =  3  Impact  Points  for  L/D  =  2.25 
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APPENDIX  D.  IMPACT  DATA  TABLES 


This  appendix  contains  all  of  the  impact  point  data  used  in  the  multiple  linear 
regression  calculations,  impact  angle  histograms  and  impact  point  scatter  plots.  The  last 
recorded  point  was  considered  to  be  the  point  at  which  the  mine  moved  outside  of  the 
background  grid's  boundaries.  The  first  four  columns  are  the  input  parameters,  t 
represents  time  and  the  last  seven  columns  are  the  output  parameters.  The  data  has  been 
organized  by  drop  angle,  decreasing  mine  length  and  decreasing  COM  position.  All  of 
the  data  contained  has  been  made  non-dimensional  by  the  conversions  listed  earlier. 
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